Role of Many-particle excitations in Coulomb Blockaded Transport 
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We discuss the role of electron-electron and electron-phonon correlations in current flow in the 
Coulomb Blockade regime, focusing specifically on nontrivial signatures arising from the break-down 
of mean-field theory. By solving transport equations directly in Fock space, we show that electron- 
electron interactions manifest as gateable excitations experimentally observed in the current-voltage 
characteristic. While these excitations might merge into an incoherent sum that allows occasional 
simplifications, a clear separation of excitations into slow 'traps' and fast 'channels' can lead to 
further novelties such as negative differential resistance, hysteresis and random telegraph signals. 
Analogous novelties for electron-phonon correlation include the breakdown of commonly anticipated 
Stokes-antiStokes intensities, and an anomalous decrease in phonon population upon heating due 
to reabsorption of emitted phonons. 

PACS numbers: 



I. INTRODUCTION 



The experimental study of electron flow through 
nanostructures has been a dynamic field of activity, with 
an eye on extending and complementing present day tran- 
sistor technologies, as well as generating entirely new 
applications. Nanoscale electronic transport spans a 
broad range of natural and artificially fabricated nano- 
structures, from carbon nanotubes, graphene nanorib- 
bons and silicon nanowires [H to spintronics and organic 
molecular electronics [!, Hj]. In particular, there has been 
enormous interest in quantum dot structures for explor- 
ing novel transport phenomena and device applications 
beyond the transistor switching paradigm, such as the ex- 
ploration of double quantum dot structures 0] for spin- 
based qubit manipulation and detection Electron 
transport through natural and artificial molecules forms 
a key research topic, especially for switching, sensing Q 
and quantum computation based applications 0]. 

In typical transport simulations, for instance in the 
widely implemented non-equilibrium Green's function 
(NEGF) formalism, it is common to include the effect 
of electron-electron or electron- vibronic interactions ap- 
proximately through an effective one-electron potential 
or self-consistent field (SCF) that needs to be computed 
sclf-consistcntly. There are many exceptions, however, 
where such an approximation may break down, especially 
when interaction energies dominate other energy scales 
of interest such as the level broadening and the device 
temperature. One such regime, well known as Coulomb 
Blockade (CB) Q, occurs when the device or channel 
capacitance is low enough that an active electron inside 
the channel can prevent a subsequent one from enter- 
ing. Such a single particle quantization of charge trans- 
fer is frequently observed in chemical reactions ||, but 



is a relative newcomer in electronic transport measure- 
ments [|| . The sequential addition of electrons in integer 
amounts disallows mean-field treatments which tend to 
smear out charges and interactions by treating all elec- 
trons on an equal footing. In contrast, solutions involve 
products of electronic occupancy and atomic displace- 
ment operators including an exclusion principle term 
that requires keeping track of every possible electronic 
or vibrational configuration through the employment of 
the many-particle Hilbcrt (Fock) space. Under these 
strong correlation conditions, energy levels must be cal- 
culated not through a simple band theory or an effec- 
tive potential, but as differences between total energies 
of the neutral and the cationic/anionic/vibronically ex- 
cited species. This is extremely difficult since it requires 
enumeration of all many-particle configurations (2^ of 
them, for N basis sets involving electronic and phononic 
coordinates!). However, such a complexity is necessary, 
as exclusion in Fock space creates a rich spectrum of ex- 
citations as well as universal scaling rules for the current 
plateau heights that are hard to cap ture a-priori using a 
modified one-electron potential [l0(. The alternate Fock 
space viewpoint (in general, a many-body density matrix 
theory) has been somewhat restricted to describe quan- 
tum dot transport [i, 11 1 and has been relatively unex- 
ploited in molecular electronics. The focus of this article 
is to illustrate the Fock space view-point of transport, 
and its experimental ramifications as far as observable 
signatures of many-particle excitations go. 

Our recent work in the area of molecular transport 
0, 14, EH, as well as compelling recent experi- 
ments triggered towards spin-based quantum computa- 
tion [B|, 0, 17, 18, [l9||, both argue for increased ac- 
tivity in this area. This article focusses on the role of 
both electron-electron and electron phonon correlations 
in non-equilibrium transport. The paper is organized into 



three broad sections. In the first section we introduce the 
Fock space viewpoint. The second section runs through 
the formalism of Fock-space transport. Section three dis- 
cusses Coulomb Blockade signatures created by electron- 
electron interactions. We elaborate on the non-trivial 
role of electronic excitations [l3[ in the interpretation of 
frequently observed I- V characteristics [2l|, [23] • It is fur- 
ther shown that electronic excitations can also result in 
intrinsic asymmetries within the channel that can provide 
an elegant approach to understanding negative differen- 
tial resistance (NDR), hysteresis effects [15|, and random 
telegraph noise j24[. The fourth section discusses the 
additional subtleties imposed by strong electron-phonon 
intereactions. We show that the inclusion of Fock space 
excitations within the electron-phonon manifold not only 
explains anomalous scaling of phonon conductance side 
bands [25|, [26| , but also predicts an anamolous tempera- 
ture distribution of the phonon population. 



II. THEORETICAL BACKGROUND 
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The Fock space approach to Coulomb Blockaded trans- 
port was proposed originally by Beenakker ll| in order 
to explain the CB conductance peak spacings and heights 
observed in semiconducting quantum dots. Figure [Hex- 
plains the difference between the one-particle and many- 
particle Fock space pictures. A set of TV single parti- 
cle energy levels generates 2 N Fock states correspond- 
ing to emptying or filling each of these level with one 
electron. In the one-electron picture, transport involves 
the addition and removal of electrons between a set of 
channel levels and two macroscopic electrodes (Fig.[TJa)). 
The levels themselves are computed by solving the one- 
electron Schrodinger equation, including electronic in- 
teractions approximately through a mean-field potential 
that modifies these levels dynamically. In the Fock space 
approach, however, the addition and removal of electrons 
leads to a transition between two entirely different multi- 
electron configurations, in this case, configurations that 
differ by a single electron (Fig. HJb)). In other words, 
electronic transport processes show up as vertical tran- 
sitions between total energies of electronic Fock states 
differing by a single electron. The situation changes a 
bit when transport involves the emission or absorption of 
phonons. In the one-particle picture, an electron jumps 
between two one-electron levels differing by the phonon 
energy (Fig. [He)), the phonon system itself driven to- 
wards equilibrium through a separate coupling to a ther- 
mal bath. In the Fock space approach Fig.QJd)), phonon- 
assisted transport shows up as horizontal transitions be- 
tween multiple copies of the electronic Fock space that 
each correspond to a different phonon number. 

In the following, we will progressively build complexity 
into our Fock space approach and attempt to identify the 
experimental ramifications of doing so. 



FIG. 1: Introducing Fock Space transport, a) A schematic of 
electronic transport in the one-particle picture. The channel 
is coupled to two electrodes which add and remove electrons 
to and from the channel, resulting in current flow. Given a set 
of N basis functions the transport problem computationally 
scales as N x N. b) In the Fock space picture such addition 
and removal processes can be understood in terms of transi- 
tions between states that differ by a single electron number. 
In such a scheme one needs to keep track of all possible con- 
figurations in the channel, thus resulting in an exponential 
2 N x 2 N scaling of the transport problem, c) In case of strong 
electron-phonon interactions, the phonon system is coupled 
to a phonon bath which is maintained at equilibrium, d) The 
inclusion of phonons calls for using copies of the electronic 
Fock space, each corresponding to a phonon number. Elec- 
tronic transport results in vertical transitions while phonon 
coupling results in horizontal ones. 



A. Single level system, SCF analysis 

Let us start with the smallest interacting system, 
namely, a real or artifical molecule with a single energy 
level capable of accomodating two spins (Fig. [5^) . The 
onsite energy of the level is eo and the Coulomb charg- 
ing energy U . The level is coupled to contacts which are 
held separately at thermal equilibrium at their respec- 
tive bias-separated electrochemical potentials ^l,r (L: 
left, R: right). A starting point for our analysis is the 
Hubbard Hamiltonian for the molecule 



H = eh + Uh^hi 



(1) 



where the operators h^i have eigenvalues and 1, while 
h = ?i| +h^. Exact diagonalizing this Hamiltonian leads 
to a Fock space consisting of four many-electron states 
(Figure [21b)), an empty zero electron state 1 00) with 
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energy 0, two one electron states 1 01) and 1 10) corre- 
sponding to an up or a down spin with energy e, and a 
doubly occupied up-down spin electron state 1 1 1 ) with 
energy 2e + U . Equilibrium occupancies of these many- 
electron states are given by the Boltzmann distribution 
P N = e -(EN-p.N)/k B T / /q j where kijT ig the thermal en- 
ergy, fi = Ep is the equi librium contact Fermi energy or 
electrochemical potential, and Q = 2 JV e" (£jY "'' jV)/ ' :fl ' r 
is the grand partition function. The average electron oc- 
cupancy is then given by (N) = J2n NPn- 

One can bypass the many-electron Fock space treat- 
ment by employing a suitable self-consistent (SCF) po- 
tential acting in the one-electron subspace, modifying the 
energies accordingly. In the spin restricted approach that 
treats both spins equally, 



Urscf = (OH/dN). 



(2) 



where (. . .) denotes a quantum mechanical average. The 
interacting term can be written as 

Hint = Uh^hi 

= U /2^2n a n s 

a 

= (U/2)Nj2^-(U/2)Y,nl 

= UN(N -l)/2, (3) 

where we have used the fact that N = n aj and = 
her, since n can only take values of zero or one. The spin 
restricted SCF potential is then given by 



U 



RSCF 



(dH mt /dN) = U((N) - 1/2) (4) 



For a given electrochemical potential, one guesses the 
value of the average occupancy (N) , uses it to calculate 
the SCF potential, and then calculates in turn the mean- 
field occupancy (N) of the level e = e + Uscf using the 
Fermi-Dirac distribution /(e) = 1/[1 + e( e- '*" feflT ], pro- 
ceeding along this line until self-consistent convergence. 

It is easy to see that the equilibrium occupancy N — fi 
[3] with respect to chemical potential (we drop the angu- 
lar term indicating average here) should be qualitatively 
different between the SCF and many-body results. In the 
former, the electron occupancy is a fractional amount, 
adiabatically changing from zero to two. In the many- 
body result, however, one does not simply multiply the 
results for one electron by two, but the electron occu- 
pancy changes abruptly between zero to one, followed 
by a plateau of width U over which the electrons are 
blockaded by the Coulomb interaction, after which the 
electron number reaches two abruptly. 

One could capture this Blockaded effect using an un- 
restricted self-consistent potential (USCF) by dictating 



that the up and down spins do not feel potentials due 
to themselves. By eliminating this self-interacting, the 
unrestricted potential for a particular spin is then given 
by 



Uuscf = {dH int /dh a ) 
= Un & 
= U{N~n a ) 



(5) 



with n = (n) and a represents the spin opposite to 
(7. This spin-dependent unrestricted potential eliminates 
the self-interaction of the level to which charge is being 
added. A self-consistent solution of the occupancy yields 
anN—fx [l3[ plot very similar to the exact result, showing 
that an unrestricted calculation can capture equilibrium 
Coulomb Blockade effects. 
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FIG. 2: Fock Space transport through a singly degenerate 
energy level, a) The singly degenerate energy level can carry 
two electrons, one of each type, b) The Fock space thus com- 
prises of 4 states. Transitions between these result in c) two 
transport channels e and e + U. Successive access of these 
transport channels results in two current jumps at different 
bias points separated by a plateau. A NEGF-SCF calculation 
is shown dotted for comparision. 



B. Single level system, Fock space transport 

Non-equilibrium turns out to be hard to mimic with 
any SCF theory, even with considerable latitude in our 
choice of the SCF potential. Let us assume the contact 
injection rates are given by 71,2/^, ignoring level broad- 
ening for the moment. Here 1,2 also refers to the left (L) 
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and right (R) contact. Both conventions are used in this 
paper to be consistent with NEGF literature. One can 
write down a master equation for a transition between 
the many-electron levels driven by the contacts 

Y^l'-.l'. (6) 

3 3 

where i,j represent the many-body states. The master 
equation, intuitively quite transparent from Figure db), 
can be formally derived by decoupling the contact and 
molecular density matrix equations in the steady-state 
limit and then invoking a Markov approximation that 
ignores memory effects such as energy-dependences in the 
contact broadening. For our simple example of the dot 
with two spin levels, the transition rates between the four 
Fock states, numbered as {|00), |01), |10), |11)}, are given 
by 

-Roo^oi = #oo->io = (71/1 + 72/2)/^- 
= #io->oo = (71/1 +72/2) /S 
Ro^n = i?io^ii = (7i/'i+72/' 2 )A 
flii^oi = i?ii^io = (7i/"'i+7 2 /"' 2 )A (7) 

where f 1>2 = 1/(1 + e ^~^/ kBT ], fa = 1 - fa, 
fa = 1/[1 + e («+^)/*» r ] 1 and f' 1<2 = 1 - fa. In 
short, electron addition processes are governed by prob- 
ability of occupancy / at the corresponding transition 
energies e or e + U by each contact electrochemical po- 
tential, while electron removal processes are governed by 
the probability of vacancy 1 — /. Thus the two current 
onset points occur at the bias situations shown schemat- 
ically in Figure [He). 

At steady-state, it is straightforward to solve these 
equations (only three of which are independent), along 
with the normalization J2i -F» = 1- The probabilities are 
then used to calculate the current injected by one contact 
(say the left one, "L') as 

I L = Y^(±e/h)[- ^ + E R % P i\ ( g ) 

i 3 3 

where the rates R L are obtained by only considering the 
individual left contact contributions to the correspond- 
ing rate, for example, -Rq ^io = 7i/l/& The ± signs 
correspond to addition/removal of electrons by the left 
contact. The resulting FV characteristic is shown in 
Pig-lid). 

In an SCF treatment, the current, shown dotted in Fig- 
ure H{d) is obtained by solving the rate equations in the 
one-electron subspace. The electron occupancy is given 
by N = (71/1 + 72/2)/ (71 +72), where the Fermi func- 
tions are evaluated at the energy e = e + Uscf- The 
SCF potential in turn depends on N as described ear- 
lier, so the calculation is done self-consistcntly. The con- 
verged energy is then used to calculate the current as 
/ = (2e/fi)7i7 2 /(7i +7 2 )[/i - /z]. 



The RSCF model that treats spins equally tends to 
give an adiabatically increasing current that reaches its 
maximum contact-dominated value 2e/h x 7172/(71+72) 
when the contact electrochemical potential fully crosses 
the level. It is important to note that charging alone 
can smear out this current, leading to a low conductance 
value spread out over a wide voltage range comparable to 
U. The RSCF potential U(N — 1/2) causes a continuous 
shift in levels with charge addition, which is accomplished 
in fractions. The unrestricted approach USCF gives an 
intervening Coulomb Blockade plateau of width U that 
separates the first spin addition (or removal) event from 
the second. The intervening 'open shell' plateau is at 
half the maximum value for complete level filling, which 
is understandable because the two spins are treated on 
an equal footing chemically, and therefore carry equal 
current. 

Compared with the SCF results above, the exact so- 
lution of the many-body rate equations reveals an inter- 
esting surprise that is actually quite illuminating. Of all 
many-electron configurations, the 1-electron states (and 
only those) are doubly degenerate, giving us a normal- 
ization condition that differs from a simpler version that 
ignores spins and simply multiplies all results by two. 
As a consequence of this sum-rule, which takes Pauli 
exclusion into account (preventing double up or down 
spin states, for example), the exact value of the open- 
shell current plateau depends on the Fermi functions 
as 7172/(71(1 + /1) + 72(1 + h)] in the large charging 
(U — » 00) limit. For sharp levels at positive bias (/1 = 1, 
f-2 = 0), this reaches 7172/(271 + 72]. For the strongly 
non- equilibrium situation corresponding to equal resis- 
tive couplings (71 = 72); the Coulomb plateau carries 
two-third of the maximum closed-shell current, in con- 
trast with the USCF result that gives a factor of half. 
This implies an interesting history dependence, in that 
the first spin added to the empty level carries more cur- 
rent than the second! If we keep track of the entire many- 
electron configuration space, it is easy to see that this 
counter- intuitive result arises because there are two ways 
of adding the first spin and only one way of adding the 
second (the other channel eliminated by exclusion) . This 
subtlety is completely washed away when we choose to 
work in a reduced N x N (or 2N x 2N for unrestricted) 
subspace instead of the full 2 2N x 2 2N configuration space. 

The SCF potential (dH/dn) was calculated by writ- 
ing the electron operators n = (h) + 5h, expanding 
the Coulomb term Uh^h^ and dropping the correla- 
tion terms 5n-\5fii completely, i.e., the Hartrce-Fock ap- 
proximation. One could include parts of the correla- 
tion term phenomenologically, by dictating that hihj « 
(1 — gij)(ni){rij), with <jy representing the exchange- 
correlation hole. This is in the spirit of Kohn-Sham 
theory, where the potential is calculated by various ap- 
proximate means. However, the effect of g is simply to 
renormalize the charging energy U that it adjoins, influ- 
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cncing at best the width, but not the height of the current 
plateaus. Thus, even for the simplest quantum dot, un- 
restricted potentials in the one-electron subspace cannot 
capture the nonequilibrium properties correctly. 

It is further illuminating when a comparison between 
the USCF and exact model is performed. Recall that 
in the USCF approach, we introduce the aspect of "self 
interaction correction" as described in Eq. [5] Here, 
the presence of an electron of a particular spin adds a 
coulomb cost to the other, but not on itself. Fig. [3Jb) 
shows the discrepancy between the USCF and the many- 
body transport calculations for a dot with two spin levels, 
coupled equally to two contacts (driving it far from equi- 
librium). As is clear from the results, the discrepancy is 
in the onset voltages, widths and heights of the various 
plateaus. While the plateau widths could be adjusted to 
fit the exact results by renormalizing the charging ener- 
gies parametrically to account for correlation effects, the 
heights are independent of these values, and depend only 
on universal factors arising from Pauli exclusion, and can- 
not be fixed in such a straightforward way, or by the usual 
DFT approach of progressively improving on correlations 
in the electronic structure. The issue is further under- 
scored by the fact that in the strongly asymmetric limit 
(71 = IOO72), where the system is essentially driven into 
equilibrium with the left contact, the agreement between 
USCF and many-body results is substantially improved. 
It is worth clarifying at this time though that in the mul- 
tilevel generalization, this correspondence becomes a lot 
harder to establish even near equilibrium, since even the 
number of current plateaus generated by a 2 N x 2 N Fock 
space approach differs substantially from its 2N x 2N 
USCF counterpart (unless the broadening functions bear 
enough poles through their energy-dependences to pre- 
cisely account for those missing conductance peaks). 

The situation gets simpler if we have asymmetric con- 
tacts 71 3> 72, so that the system approaches equilibrium 
with the left contact. As Fig. [3] (d) shows, the USCF 
agrees with the many-body limit in this asymmetric cou- 
pling case, as both approaches are dealing with a near- 
equilibrium problem. For positive bias on the weaker 
contact, the stronger contact keeps the level filled (which 
it can do in two ways, adding an up or a down spin, 
assuming the level was empty to begin with). For oppo- 
site bias, the stronger contact empties this level, which 
can now be done in only one way (up OR down depend- 
ing on what occupied the level). The competition be- 
tween 'shell tunneling' and 'shell filling' makes the I-V 
strongly asymmetric, with the first plateau half the sec- 
ond for positive bias on the weaker contact, and merging 
with the second for opposite bias. This ratio of one to 
two, observed experimentally (27j , arises in a straightfor- 
ward way from our analyses since the ratio of the first 
and second plateau currents is given for positive bias 
by (71 + 72)/(27i + 72) ~ 1/2 for 71 > 72, and by 
(71 + 72)/(7i + 272) ~ 1 for negative bias. The asym- 



metry arises from the difference in the number of spin 
addition and removal channels for positive and negative 
bias, and leads to an asymmetry in the current levels. 
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FIG. 3: Comparison between USCF and exact transport re- 
sults for (a) extreme non-equilibrium (71 =72), b) in which 
the I-V's show significant discrepancy between USCF (dot- 
ted red) and exact (bold black) result, c) The case of strong 
asymmetry implies a near equilibrium with the first contact 
(71 = IOO72). d) The discrepancies in plateau onset, width 
and height are resolved in this equilibrium limit, but are hard 
to resolve once we go far from equilibrium. 

We will next show how this model is extended for a 
larger molecule, and how additional physics due to cor- 
relations and excitations start to arise. 



C. General Approach for multilevel systems 

In case of a larger molecule, one begins with the model 
Hamiltonian in second quantized notation: 



H = e a h a 



a 
a, a 



W) c i c /3 



(9) 



where h a = c* a c a , a, (3 correspond to the orbital indices of 
the orbitals for various sites on the molecule, and a, a rep- 
resent a particular spin and its reverse. Exact diagonal- 
izing this Hamiltonian yields a large spectrum of closely 
spaced excitations in every charged molecular configura- 
tion. Using the equation of motion of the density matrix 
of the composite molecule and leads and assuming no 
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molecule-lead correlations, one can derive [28|, [29| a sim- 
ple master equation for the density-matrix of the system. 
Ignoring off- diag onal coherences, we are left with a mas- 
ter equation [28| in terms of the occupation probabilities 
Pf of each N electron many-body state \N, i) with total 
energy . The master equation then involves transi- 
tion rates R(N,i)^>(N±i,j) between states differing by a 
single electron, leading to a set of independent equations 
defined by the size of the Fock space [ll | 

\p pN p pW±ll 



dt 



E 

N.j 



l(JV±l,j)-(JV,i)-fj 



along with the normalization equation J^. N P i 



N 



(10) 

= 1. 

For weakly coupled dispcrsionless contacts, parameter- 
ized using bare-electron tunneling rates j a , (a: left /right 
contact), we define rate constants 

r$; = la \{N,i\ct\N-i,j)\* 

?iia = j a \(N,i\c a \N + l,j}\ 2 , (11) 



c* a , c a are the creation/annihilation operators for an elec- 
tron on the molecular end atom coupled with the corre- 
sponding electrode. The transition rates are given by 



R 



(N,i)->(N-l,j) 



R 



(N-l,j)->(N,i) 



E 

a=L,R 

E 

a=L,R 



■)] 



ij 



r f(e 



Nr 
ij 



Mo 



(12) 



for the removal levels (N,i — » N — 1, j), and replacing 
(r — ► a, f — ► 1 — /) for the addition levels (N, i — > N + 
l,j). n a are the contact electrochemical potentials, / 
is the corresponding Fermi function, with single particle 
removal and addition transport channels e 



Nr 



E 



N 



£f-\and e£ a 



E 



N+l 



E™ . Finally, the steady-state 



solution to Eq. (TIT)]) is used to get the left terminal current 



I = ±- 



R 



(N,i)- 



pN pL 



■(N,i)Pj 



N±l 



(13) 

where states corresponding to a removal of electrons by 
the left electrode involve a negative sign. We will assume 
a break-junction configuration with equal electrostatic 
coupling with the leads, ^l,r = E F =p eVd/2. 

While the above equations include spectral details from 
the multiple excitations, a considerable simplification 
arises if we can incoherently sum over many of these ex- 
citations, leading to the 'orthodox model', where 

1 = E [ R N^N±l - R L N^N T l] P N (14) 
n N 

The transition energies can be obtained in terms of a sim- 
ple RC circuit, while the transition rates also compactify 
once we integrate them over the relevant energies, taking 
exclusion factors into account. 

The implementation of the above sets of equations now 
sets the stage for further discussions. 



III. COULOMB BLOCKADE: THE CASE OF 
ELECTRONIC TRANSPORT CHANNELS 



In this section we describe how the excitation spec- 
tra of molecules may be probed using Coulomb Block- 
ade transport. While typical molecular I-Vs look rela- 
tively featureless, what is not often appreciated is that 
in the Coulomb Blockade limit (realized by engineering 
weak contacts or non-conductive backbones), electronic 
excitations do give rise to prominent features in the I-V 



211 . [22| . We will also show how such an excitation spec- 
tra can generate intrinsic I-V asymmetries, which in turn 
can manifest itself as NDR effects. 

Exact-diagonalization of the molecular Hamiltonian 
Eq. [5] generates a body of excitation spectra for every 
charged configuration. Our typical starting point is the 
equilibrium molecular quantum dot configuration in its 
ground state. Addition or removal of an electron takes 
the system to the a new ground state corresponding to 
the singly charged cation or anion, marking the onset 
of conduction. In addition to the ground state however, 
each charged species bears a quasi continuous excitation 
spectrum separated from its ground-state energy by a 
gap that is determined by the energetics of the molecular 
system. Once the first excitation is accessed, the quasi- 
continous excitation spectrum can be easily probed, as 
shown in Figures SJa), (b) and (c). It is worth empha- 
sizing that the excitations are energetically quite close as 
they differ in the rearrangement (but not in the number) 
of charges. At high bias with small broadening and in 
the absence of phonons, the purely electronic excitations 
do not have adequate time to relax to the ground state, 
accounting for their visibility in the current spectrum. 

The simplest and most prominent impact of Coulomb 
Blockade on the TVs of short molecular wires is a clear 
suppression of zero-bias conductance, often seen exper- 
imentally [3, However, integer charge transfer can 
also occur between various electronic excitations of the 
neutral and singly charged species at marginal correla- 
tion costs [U [33 ■ The above fact leads to fine structure 
in the plateau regions 2(| 21 , 23|, HI] , specifically, a quasi- 
linear regime resulting from very closely spaced transport 
channels (e^) via excitations. The crucial step is the ac- 
cess of the first excited state via channel e^ r , following 
which transport channels involving higher excitations are 
accessible in a very small bias window. 

The sequence of access of transport channels upon bias, 
enumerated in the state transition diagrams shown in 
Figs. |Ua),(b) and (c) determines the shape of the I-V 
characteristic. When the Fermi energy Ep lies closer 
to the threshold transport channel €qq corresponding to 
charge transfer between two ground states, it takes an 
additional positive drain bias for the source to access the 
first excited state of the neutral system via the transi- 
tion €i a r . Under this condition, the I-V shows a sharp 
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FIG. 4: Electronic excitations in molecular conduction, a) 
Schematic shows three Fock states, the N — n equilibrium 
ground state configuration and two states in the addition 
N — n + 1 spectrum, b) Each charge configuration consists 
of a ground state and a set of closely spaced excitations. The 
access of such closely spaced excitations between the neutral 
and char ged species gives rise to numerous c) transport chan- 
nels e trk hj. The access of such transport channels gives rise 
to a quasi ohmic behavior in the I-V characteristics that are 
commonly noted in various experiments [2(], H3| . 



rise followed by a plateau (Fi gure Efd), shown in blue), 
as seen in various experiments [321 ]. However, when trans- 
port channels that involve low lying excitations such as 



ef r are closer to the Fermi energy Ep than Cqq , the ex- 
citations get populated by the left contact immediately 
when the right contact intersects the threshold channel 
e^f, allowing for a simultaneous population of both the 
ground and first excited states via €qJ and e^ r at thresh- 
old. Under these conditions the I-V shows a sharp onset 
followed immediately by a quasilinear regime (Fig. SJd) 
shown in black) with no intervening plateaus, as observed 
frequently in FVs of molecules weakly coupled with a 
backbone USED,!!!. 



A. Extrinsic Asymmmetries -the case of rectification 

The direct role of excitations in conduction becomes 
particularly striking under asymmetric coupling (7^ >> 
7ft) with contacts (2ll.[35|. Due to this extrinsic asymme- 
try, current magnitudes are dictated by the weaker con- 
tact. This asymmetry directly affects the forward and 
reverse bias characterestics, leading to current rectifica- 



FIG. 5: Current asymmetry from extrinsic "contact" asym- 
metry, a) In the case when the left electrode is more strongly 
coupled to the channel than the right one 7_l >> 7h, cur- 
rents are limited by removal (addition) rates under forward 
(reverse) bias conditions, b) The number of removal chan- 
nels can be significantly different from the number of c) ad- 
dition channels. This gives rise to the asymmetry between d) 
forward and reverse bias situations. Furthermore the appli- 
cation of gate potential modulates the number of transport 
channels available at threshold. This feature results from the 
Fock space excitation spectra, causing not only a change in 
threshold voltage with applied gate bias, but also a discern- 
able threshold current modulation as noted here, and in var- 
ious experiments [2lj . 



tion. This rectification is caused by the inherent asym- 
metry between addition and removal processes, each of 
which is the rate limiting depending on the bias direction, 
as shown in Figs.^a), (b) and (c). 

In contrast to the SCF regime where unequal charg- 
ing drags out a same level current over different volt- 
age widths (36[, in the CB regime we encounter clear 
intermediate current steps from open shells, with current 
heights that are themselves are asymmetric at threshold 
(Fig. [5](d)) . This asymmetry arises due to the difference 
in the number of pathways for removing or adding a spin, 
taking in particular into account the possible excitation 
channels between the neutral and singly charged species 
(Figures [5]Jb) and (c)). The number of such accessible 
excitations at threshold can be altered with an external 
gate bias, leading to a prominent gate modulation of the 
threshold current heights, over and above the modula- 
tion of the onset voltages and the conductance gap plj 
(Figure EJd)). Furthermore, it is easy to show that the 



8 



asymmetry will flip between gate voltages on cither side 
of the charge degeneracy point, as is also observed exper- 
imentally [35j] . 

It is worth emphasizing that the sophistication arose 
specifically due to the presence of separate identifiable 
contributions from the open shells, which are normally 
overwhelmed by broadening if we were away from the CB 
limit. However, the specific identification of these exci- 
tations is not crucial to the qualitative shape of the I-Vs, 
as long as one is not looking too closely at the individ- 
ual spectral features. A simpler orthodox model would 
then suffice to reproduce these broad features described 
above, such as the transition between steps and slopes in 
the I-Vs, the flipping of asymmetry and the gate modula- 



tion of the current levels [14|. However, there are notable 



exceptions, such as intrinsic asymmetries, where a clear 
separation needs to be made between certain classes of 
excitations with longer lifetimes ('traps') and the regu- 
lar excitations with shorter lifetimes ('channels') that are 
responsible for the current flow. 



B. Intrinsic Asymmetries - NDR, hysteresis and 
telegraph noise 

The physics of NDR can be explained readily using an 
USCF model that is actually quite intuitive. Consider a 
channel and a trap, the channel being strongly coupled to 
the contacts and the trap weakly coupled. Accessing the 
channel creates a resonant onset of current. Accessing the 
trap subsequently would keep the conduction unaffected, 
as the trap is non-conducting; however, once we throw in 
the strong Coulomb repulsion arising from charging up 
the trap, it is easy to sec that this repulsion can expel the 
conducting channel out of the conducting bias window, 
leading to an NDR. Self-interaction correction is crucial 
to this picture, as the charging should repel the channel 
level, but not the trap level itself. Further sophistications 
can arise by considering the lifetime of the trap. If the 
trap does not release its charge during the measurement 
time, then a reverse scan would keep the channel blocked 
and lead to a hysteresis. Such a hysteresis is scan-rate 
dependent, as the rate determines the degree to which 
the trap releases its captive charge. The low current state 
can be reversed by going to large negative bias to expel 
this charge and restore initial conditions. 

The realization of the NDR involves two conditions: 
(a) that the charging of the trap is large enough to expel 
much of the channel from the bias window, and (b) that 
the resulting current carried by the trap plus the residual 
tail of the channel sitting in the bias window exceeds the 
current carried originally by the channel alone. The first 
amounts to the onset condition for an NDR, while the 
second describes, in some sense, the effectiveness of the 
NDR (in terms of an inequality). With a little bit of 
effort, one can extract a range of parameters that satisfy 



both these conditions. We will instead show that even in 
the Fock space picture (which treats the channel plus trap 
as a composite system), one can identify these two same 
conditions: the NDR starts when one encounters what 
we refer to as a 'dark' or 'blocked' state, and subject to 
this condition, the effectiveness of the NDR amounts to 
an inequality that we will now discuss. 

A condition for NDR can be derived much more gener- 
ally in terms of three Fock space states \A), \B), and \C) 
with energies Ea < Eb < Ec respectively (FigH^a)), 
representing three accessible states within the bias range 
of interest. For instance, \A), \B) could be the ground 
states of the N = no and N = no + 1 electron systems, 
while \C), the first excited state of the N = no + 1 
electron system. Transport of electrons involves single 
charge removal or addition between states \B), \C) and 
|^4) that differ by an electron, via addition and removal 
rates Ra^b c oc — - — . Such an electron exchange is 

^ <^ TAB,C ° 

initiated when reservoir levels are in resonance with the 
single electron transport channels 6ba = Eb — Ea and 
cca = Ec — Ea respectively. The I-V characteristic of 
this three state system, shown schematically in Fig[BJb), 
shows two plateaus with current magnitudes I a and I& re- 
spectively. Current collapse or NDR occurs when I a > I b . 

In the above system, NDR occurs when under specific 
conditions, the state \C) can be a blocking or dark state, 
for which electron addition is feasible while removal is 
rate limiting. This can happen when there arc intrin- 
sic asymmetries within the transport problem, as in the 
trap-channel dichotomy described earlier. Regardless of 
the specific origin of blocking, a simple criteria for current 
collapse can be derived [l5| in terms of the rate limiting 
process, that is the electron removal rate Rc^A from the 
dark state \C). A better intuition is provided by think- 
ing of rates as inverse lifetimes, as shown in Figure [|3c). 
If the life time of state |C), tca exceeds the sum of the 
addition and removal rates of the conducting state, it 
becomes an effective blocker. In other words, the condi- 
tion Tqa > t ab + t ba determines whether NDR occurs 
or not. The superscripts L and R represent left or right 
electrode, and under the forward bias situation add and 
remove electrons respectively. The consequence of the 
above criterion is summarized in Figure [61(d), clearly in- 
dicating that NDR only occurs when tq A > t\ b + t§a- 
A similar criteria is valid for the reverse bias direction by 
interchanging superscripts L with R. 

Further sophistications may arise due to 'transient' 
probing. The criteria obtained above applies for steady 
state, where the dark state gets occupied with certainty, 
resulting in an inevitable current blockade regimee. Let 
us consider such a blocked current state, shown dotted 
in Figure Etc). A 'disguised' NDR can also be achieved 
if the voltage sweep rates are faster than the rate deter- 
mining dark state life time. Such a sweep rate simply 
means that the dark state has not yet been occupied and 
thus a current determined by the addition and removal 
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FIG. 6: NDR effects from intrinsic asymmetry or "dark 
states", a) A generic mechanism NDR can be cast in terms 
of three device Fock-space states, and transitions | A) B) 
and | A) <->| C) between those that differ by a single elec- 
tron. The dark state | C) say has a very slow removal rate 
in comparison to its addition rate, b) Bias environments, a 
and b, correspond to current rise followed by a collapse, c) A 
general criteria for such a current collapse or NDR to occur 
with increasing voltage can be cast in terms of the dark state 
removal rates, d) NDR in the I-V characteristic is contingent 
to the criteria derived here Hal. 



times l/(T~A B + Tg A ), can still flow. This situation is 
shown in Figure [T^c), in which the onset of dark state 
can be delayed provided the sweep rate is fast enough. 
But once the dark state forms, the current is blocked 
and thus remains so. This implies that on the backward 
sweep current remains blocked, resulting in a sweep rate- 
dependent hysteresis. Such a behavior has been noted 
by Kiehl et.al. j3?J in a molecular system, which was at- 
tributed to a slow charge trapping process, which again 
fits into the dark state picture. Just as noted here, Kiehl 
et.al. also observe that the NDR effect diminishes with 
decreasing sweep rate and ultimately vanishes at steady 
state, implying the absence of a real hysteresis at steady 
state. This kind of hysteresis is just an artifact of sweep 
rates being faster than the rate determining processes. 

It must be mentioned that certain hysteretic processes 
can occur even at steady state due to the true bistable 
nature of the system. A classic example is the charg- 
ing induced hysteresis in resonant tunneling diode (RTD) 
structures [381 ] . Other interesting examples include bista- 
bility caused by hyperfine interactions [40( causing a hys- 
teresis with applied magnetic field, and optical bistability 



FIG. 7: Hysteresis effects due to sweep rate, a) Under steady 
state conditions the dark state behaves as an effective blocker, 
b) A sweep rate faster than the rate determining time scale 
can put the transport process in a transient state, c) Un- 
der such transient conditions sweep rate can induce disguised 
NDRs and often hysteresis, both depending on sweep rate 

51. 



created by superlattice dielectrics that show strong tera- 
hertz nonlincarity due to Bloch oscillation and dynamic 
localization of electrons 4l|, |42| • 

While we discussed the interesting case of dark states 
earlier, touching upon relevant consequences such as 
NDR and hysteresis, we avoided particular examples. Re- 
cently, we applied our dark state model [l5| to explain 
the NDR observed due to subtle spin correlation effects 
in double quantum dots 1391] . Indeed this NDR results 
in the spin blockade regime [391 ] which currently forms a 
key concept in the area of single spin manipulation and 
control @, [3, [H El, 

Finally, the transient probing with a high resolution 
probe would allow us to explore the approach to reso- 
nance with the trap states. As we discuss in 24|, the 
stochastic blocking and unblocking of the channel by 
the occupation/deoccupation of the traps near resonance 
generates a flicker in the output current known as ran- 
dom telegraph noise. The ratio of capture and emission 
times is given by a Boltzmann factor whose argument 
depends on the trap energy location. Shifting this trap 
with a combination of gate and drain voltages leads to an 
associated scaling of the capture to emission time ratio, 
allowing one to infer the spatial and spectral location of 
the trap states. This provides a powerful 'barcode' for 
characterizing single-molecular defects. 
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Our analyses over the last few sections focused on 
electron-electron correlations and their treatment in a 
Fock space approach that captures the relevant physics 
(if at the expense of computational simplicity). In the 
next section, we will discuss how this approach can be 
extended to incorporate electron-phonon correlations in 
transport. 



IV. ELECTRON-PHONON INTERACTIONS: 

So far, our primary concern was Coulomb Blockade 
and electronic excitations. Here, we show that coupling 
Coulomb Blockade with phonon-assisted tunneling re- 
sults in non-trivial physics. This section is primarily mo- 
tivated by a recent series of experiments performed on 
suspended carbon nanotube quantum dots [HI]. In such 
a suspended system, the phonons are driven far out of 
equilibrium and couple strongly with the electronic sys- 
tem. Our approach then is to invoke the Fock space of 
the electron phonon system as indicated in Figures HJc) 
and (d). 



A. SCF treatment — IETS and phonon sidebands 

The SCF approach to electron-electron interactions 
has been discussed at length in our earlier papers and 
contrasted with the Fock space approach. It is worth 
quickly touching upon the SCF treatment of phonons, 
before diving into its Fock-space analogue. In presence 
of dephasing scattering events, the current at the left 
contact can be written in the NEGF formalism [43| as 



^ J dETr[Y,fA - T L G r - 



(15) 



where A is the spectral function, is the broadening 
by the left contact, E™ is the in-scattering self-energy 
from the left contact and G n is the correlation func- 
tion describing the energy-dependent occupancy of the 
levels, taking quantum interference into account. The 
influence of scattering by contacts and phonons sits in 
£ m = ££» + £™ + . The equations for the contact Ss 
are well known - we will focus here on the additional con- 
tributions from the phonon scattering processes. Within 
the self-consistent Born approximation (self-consistency 
needed to conserve current), the equations connecting 
the phonon contributions form two groups - a set of dy- 
namic equations, and a set of static equations. The non- 
equilibrium dynamics describing the filling and emptying 
of states is described by 

G n ' p (E) = G{E)^ out {E)G\E) 



(E) = D (ui) ® 



(16) 



where ne(w) = [ e h "/ k BT _ jj-i j g |- ne equilibrium Bose- 
Einstein distribution of the phonons, and Dq is its de- 
formation potential, in other words, (-Do)ij — hihj, h a 
being the electron phonon coupling constant at the real 
space point a (In general, ft is a non-diagonal matrix 
in an arbitrary basis set, and Dq is a fourth-rank ten- 
sor). ® denotes an element-by-element (as opposed to 
matrix) multiplication. The equations above are for a 
single phonon mode at frequency to, and will need to be 
integrated over a phonon density of states for a continu- 
ous distribution of phonons. 

The static equations describing the states themselves 
are given by 



T ph (E) 
G(E) 



z;ke) 

H(X ph )+i 



ph 

2 



[EI-H - S L - S fl - S p7l ]- 



(17) 



n B (w)G"' p (£TM 



where TL denotes the Hilbcrt transform. The above equa- 
tions capture the physics of phonon-assisted tunneling 
in all its various limits. In its simplest form, it con- 
tributes to dephasing that reduces the ballistic current 
in devices. Near resonance, the ±fko terms in the ar- 
guments of the G n,p matrices contributing to T, m ' out (E) 
generate phonon sidebands. The scaling of these side- 
bands depends on the deformation potential as well as 
the distribution of the phonons. We have assumed this 
to be equilibrium Bosc-Einstcin rig, although one could 
generalize it provided we have a separate evolution equa- 
tion for the phonon dynamics coupled to the electron 
transport equations that describe parametrically or oth- 
erwise how the phonons are driven away from equilibrium 
by the electronic subsystem, and how the other terms in 
the phonon evolution equation (typically involving cou- 
pling with a thermal reservoir) try to bring this back 
to equilibrium. In the next section, we will, in fact, in- 
clude these 'hot' phonon equations explicitly; instead of 
separate coupled equations, however, we will treat the 
electron-phonon as a composite system, and use the cou- 
plings between the electrons, phonons, and with the con- 
tacts and the bath as processes driving the evolution of 
the composite system. 

The NEGF formalism also allows us to get away from 
this phonon-assisted tunneling limit to the off-resonant 
limit. When the electronic levels lie far from resonance, 
their sidebands do not show up and we get instead the 
inelastic electronic tunneling spectrum (IETS) of the 
molecular system. For weak electron-phonon coupling 
(retaining leading order terms in Do), the NEGF algebra 
simplifies considerably, so that the current partitions into 
an elastic component given by the usual Landauer for- 
mula, as well as an inelastic component that explicitly in- 
volves exclusion principle terms at different emission and 
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absorption energies. For applied voltages larger than the 
phonon frequency, these terms create additional current 
transport channels through phonon emission, creating a 
slight increase in the current that only shows up as peaks 
in the second derivative with voltage, generating the fa- 
miliar IETS spectrum. The physics resides entirely in 
the inelastic current - one can use sophisticated electron 
structure methods to extract these peak positions from 
the phonon frequencies, as well as the IETS peak heights 
from the computed deformation potentials. The results 
also show additional subtleties near resonance that are 
experimentally observed. Specifically, near resonance the 
clastic current also picks up phonon signatures from the 
phonon contributions to S residing in G, generating a dip 
rather than a peak that arises from the phonon sidebands 
described above (the equations also generate a familiar 
polaronic shift of the main peak). The inclusion of the 
equilibrium ns distribution gives phonon emission peaks 
that are stronger in strength than absorption peak by 
Boltzmann ratios evaluated at the electronic tempera- 
ture. 

While encouraging agreement between computed and 
experimental off-resonant behavior, specifically, IETS 
spectra has been reported [3| , the resonant phonon side- 
bands and the scaling of their conductance peak heights 
with current shows significantly more complex behavior 
arising from strong non-equilibrium clcctron-phonon cor- 
relations, which necessitates a Fock space approach. We 
will now introduce the joint electron-phonon Fock space 
and the experimental ramifications of electron-phonon 
correlation effects captured with this treatment. 



B. Fock space treatment of phonon sidebands 
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FIG. 8: Fock space model to include electron-phonon interac- 
tions, a) The dot is electrically connected to the left (right) 
contact (with electron tunneling rates Tl,r/?i) and mechani- 
cally to the phonon bath (with a phonon escape rate of /3/h). 
The dot has electronic degrees of freedom and phonon de- 
grees of freedom huij = 1,2,3,...) with coupling Xij. 
b) Transition between different states of a dot, having one 
electronic level and one phonon mode, due to coupling to 
the contact and the phonon bath (|0,iV > and \1,N > de- 
notes the states of the system having N phonons each but 
electron and 1 electron respectively), c) Effective Coulomb 
blockaded channel in the absence of phonons. d) Inclusion 
of electron-phonon coupling adds phonon side bands to the 
existing transport channels. 



We start from a model Hamiltonian for a quantum 
dot having onsite energies e,, Coulomb interaction en- 
ergy Uw , vibronic modes at energy fkOj and electron- 
phonon coupling . The system is connected to electri- 
cal contacts with couplings Ti^ and to a thermal bath 
with coupling (3 (figure (Ha)). Electronic transport due 
to bias applied to the contacts as well as phonon emis- 
sion and absorption processes lead to transitions between 
many-body states \e l Ne , k) (k phonons and ith electronic 
level in the N e electronic subspace) of the quantum dot 
(figure Efb)). The rates of these transitions due to left 
(right) contact (TZ^ r k)->\e s p)) an( ^ * ne phonon bath 



,) are calculated by applying Fermi's 



Golden Rule [26| . When an electron is added or removed 
via standard transport processes, it can also change the 
number of phonons, taking the dot from a state |0,iV) 
to 1, N ± p) or vice versa, where p is the number of 
phonons emitted or absorbed. When the quantum dot 
absorbs or emits a phonon, the state changes from | X, N) 
to \X, N ± 1), where X = or 1. The consequence of 



phonons in transport is the addition of extra transport 
channels to the already existing Coulomb Blockade trans- 
port channels, as shown in Figures [S^c) and (d). 

The rate equations in the above Fock space ([5]) can 
be solved for a finite number of phonon emission or ab- 
sorption channels to extract the resulting current-voltage 
characteristic. Focusing on the intriguing experiments on 
Coulomb Blockaded nanotube quantum dots with promi- 
nent breathing modes 25j, one sees multiple intrigu- 



ing features: (a) the absorption and emission sideband 
heights do not scale as simply as above. This is because 
the number of phonons N p h itself changes with current, 
in addition to being driven far from equilibrium in the 
suspended sections of the tube with small escape rate 
f3. (b) The scaling of sidebands with phonon popula- 
tion differs significantly from predictions of the analogous 
Tien-Gordon theory of photon-assisted tunneling (PAT) 
[Ifil] . This discrepancy arises because unlike the PAT ex- 
periments, the phonons are not coherent, and are partly 
correlated with the nanotube electrons ]2&i 
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C. Non-equilibrium phonons and effective 

temperature 

A crucial component of the above picture is the 
strongly non-equilibrium distribution for the phonon 
population, as their generation rate by the drive cur- 
rent (determined by Tl,r and A) exceeds their extrac- 
tion rate (3. From the solution to the rate equations 
giving us the joint electron-phonon occupation proba- 
bilities, we can calculate the probability distribution of 
the phonon by summing over the electronic subspaces: 
~ Ne P\e" N ,k) ■ Defining an effective temperature 
T* corresponding to the non-equilibrium phonon occu- 
pation N p h = [e hul l kBT — 1] _1 , we find the correspond- 
ing Boltzmann distribution: Pf 1 ® = e -khw/k B T* j z of 
phonon subspace occupation probability, where Z is the 
partition function (T* is identified by fitting the higher 
energy tails between the two distributions). A compari- 
son between P£ h and P% hB at different bias for different 
decay rates of the phonons reveals that they differ con- 
siderably as the phonon decay rate (3/h decreases (fig- 
ure Efa)) at some applied bias. The deviation from a 
Bose-Einstein like shape suggests that the phonons in 
the suspended nanotubes are strongly non-equilibrium, 
so that temperature is not a well-defined quantity except 
to describe the higher energy tails (figure [9f b)). 

D. Effect of the surrounding temperature: 

anomalous phonon population 

The non-equilibrium phonon population reveals a 
striking and somewhat counter-intuitive dependence on 
the substrate temperature, the phonon population de- 
creasing with increasing temperature at certain bias volt- 
ages (Fig.fTOTa')'). This peculiar behaviour manifests itself 
as long as the surrounding temperature is smaller than 
the separation between the upper emission and lower ab- 
sorption sidebands of two subsequent Coulomb Block- 
ade peaks. Beyond this temperature the phonon occu- 
pancy increases monotonically with temperature for all 
bias values, as expected. The anomalous temperature 
dependence arises from a trade-off between phonon gen- 
eration and recombination rates in the coupled dot-lead- 
bath system. The behavior is observed specifically at 
bias values corresponding to the onset of a new phonon 
absorption channel (Figure. flffl b). bottom left). With 
increasing temperature from ~ 5 to 20K, the increas- 
ing tails of the contact Fermi functions redistribute the 
electrons from a phonon emission sideband t\ + huj of a 
lower electronic peak to an absorption sideband £2 — hw 
of a higher peak (Figure. ITOTb) bottom right). Under 
this specific conspiracy of temperature and bias values, 
the number of electrons resonant with the phonon side- 
bands decrease with increasing temperature, decreasing 
the efficiency of the phonon-assistcd transport. At higher 
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FIG. 9: The phonon distribution P£ from the Boltzmann 
dist. &)(solid line) with the effective temperature T* differs 
considerably from the phonon distribution P£ h calculated di- 
rectly from the rate equation (circles), b) The phonon dis- 
tribution P£ h approaches the Boltzmann distribution P£ hB 
(with proper effective temperature T* not shown here) at el- 
evated surrounding temperature T. 

temperatures between 20-50 K, the temperature is large 
enough to open new emission channels that eventually 
increases the phonon population at all bias voltages. 

CONCLUDING REMARKS 

For a wide variety of transport problems, a perturba- 
tive treatment of interactions coupled with a quantum 
kinetic theory (such as NEGF) does an admirable job of 
explaining and predicting experimental features and pro- 
viding important physical insights. As objects scale to- 
wards nanodimensions, however, strong confinement and 
poor coupling with the surroundings lead to increasing 
degrees of correlation, especially at lower temperatures. 
In this limit, the Fock space approach (more generally, 
the many-body density matrix approach) naturally al- 
lows us to compute transport signatures, provided we 
have a suitable means of extracting the various corre- 
lated many-body states and the transition rates among 
them. Exact diagonalization provides one option, al- 
though this quickly becomes computationally intractable. 
Partial configuration interaction (CI) schemes may prove 
to be more practical. In contrast with the NEGF-SCF 
limit where one could aim for quantitative and predictive 
accuracy with increasing amounts of chemical sophistica- 
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FIG. 10: Anamolous temperature dependence of phonon dis- 
tribution a) The number of phonons decreases at some bias 
points with increasing surrounding temperature T. (b) Bot- 
tom LeftAt lower surrounding temperature excitation of en- 
ergy ei + hui associated with resonant level ei falls inside the 
bias window between \il + UbTl and fXR — ksTu and the ex- 
citation of energy £2 — hui associated with resonant level £2 
falls outside the bias window. Bottom, right: At elevated sur- 
rounding temperature fermi functions of the contacts broad- 
ens around respective chemical potentials and the excitation 
of energy £2 — hui enters the bias window and lesser no. of 
electrons appear at energy £i + huj in the left contact. 



tion, such details are hard to build into the Fock space 
approach and need to be replaced by simpler, model 
problems, with parameters that could be benchmarked 
with more detailed models and measurements. However, 
even these simple models (with a reasonable choice of pa- 
rameters) show qualitatively new physics that is exper- 
imentally observable, ranging from gate-modulated cur- 
rent levels, gate tunable excitation spectra, scan-rate de- 
pendent NDR and hysteresis, and the breakdown of our 
common intuition based on equilibrium phonon sideband 
scaling and the classical theories of photon-assisted tun- 
neling. While these two limits (quantum wire and quan- 
tum dot) are separately well understood, at least for- 
mally, the intermediate coupling regime between the two 
becomes particularly challenging to model as there is no 
small 'fine structure' parameter that would allow a con- 
venient starting point for a perturbation expansion (e.g. 
a nonintcracting wire or a fully interacting but isolated 
quantum dot). Significant progress is needed at formal, 
computational and experimental levels in order to probe 
this regime, which bears the promise of completely novel 



physics of non-equilibrium correlations as well as possible 
applications based on the interaction between conducting 
detector elements in the quantum wire regime and non- 
conducting storage elements in the quantum dot regime. 
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